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The IceCube neutrino detector is a cubic kilometer TeV to PeV neutrino detector under con- 
struction at the geographic South Pole. The dominant population of neutrinos detected in IceCube 
is due to meson decay in cosmic-ray air showers. These atmospheric neutrinos are relatively well- 
understood and serve as a calibration and verification tool for the new detector. In 2006, the detector 
was approximately 10% completed, and we report on data acquired from the detector in this con- 
figuration. We observe an atmospheric neutrino signal consistent with expectations, demonstrating 
that the IceCube detector is capable of identifying neutrino events. In the first 137.4 days of live- 
time, 234 neutrino candidates were selected with an expectation of 211 ± 76.1(syst.) ± 14.5(stat.) 
events from atmospheric neutrinos. 

PACS numbers: 95.55. Vj, 95.85.Ry, 96.50.sf 



I. THE ICECUBE DETECTOR 

The IceCube neutrino detector is being deployed in 
the deep ice below the geographic South Pole [l(. When 
completed, the detector will consist of two components. 
The Inlce detector is a cubic kilometer of instrumented 
ice between 1.5 and 2.5 kilometers below the surface. A 
cubic-kilometer has been long noted as the required scale 
to detect astrophysical neutrino sources above the atmo- 
spheric neutrino background (see, e.g. [2] 0] 0)- T nc 
IceTop detector is a square-kilometer air-shower array at 
the surface. This analysis concerns data from the Inlce 
detector exclusively. 

The Inlce detector consists of an array of light-sensitive 
Digital Optical Modules (DOMs) |[, deployed 17 me- 
ters apart in strings of 60. Strings are arranged on a 
hexagonal grid with a spacing of 125 meters. The DOMs 
house a 10-inch Photomultiplier Tube (PMT) and elec- 
tronics to acquire, digitize and time stamp pulse wave- 
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forms from the PMT. With a waveform fit for fine tim- 
ing, the timing resolution for individual photon arrivals 
is expected to be less than 2 nanoseconds [f|. In 2006, 
the DOMs were operated in Local Coincidence (LC) with 
their neighbors, meaning that a triggered DOM's wave- 
form was only transmitted to the surface if an adjacent 
DOM on the string also triggered within ±1000 ns. The 
surface data acquisition system forms triggered detector- 
wide events if 8 or more DOMs are read out in 5 /.is. 

The detector is being deployed in stages during austral 
summers from 2004 to 2011. After the 2005-2006 season, 
the Inlce detector consisted of 9 strings, termed IC-9. 

IceCube is optimized for the detection of muon neu- 
trinos in the TeV to PeV energy range. It is sensi- 
tive to these muon neutrinos (and muon anti-neutrinos) 
by detecting Cherenkov light from the secondary muon 
produced when the neutrino interacts in or near the in- 
strumented volume. Neutrino-induced muons are sepa- 
rated from air-shower-induced muons by looking only for 
muons moving upward through the detector. Up-going 
muon events must be the product of a neutrino interac- 
tion near the detector, since neutrinos are the only known 
particles that can traverse the Earth without interacting. 



II. ATMOSPHERIC NEUTRINOS 

Neutrinos produced in cosmic-ray air showers at the 
Earth arc known as atmospheric neutrinos and form the 
chief background to potential astrophysical neutrino ob- 
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servation. The atmospheric neutrino spectrum is rela- 
tively well-understood [6] 7] and has been measured up 
to 10 5 GeV by AMANDA 0. Atmospheric neutrinos 
from the decay of charmed mesons can contribute signif- 
icantly above 10 4 GeV, depending on the model (see e.g. 
iPildll). This prompt component is not well known 
due to uncertainties in the charmed meson production, 
but with the present exposure of IC-9, this prompt com- 
ponent is negligible and it is presently neglected. 



III. RESULTS 

Data acquired from the IC-9 detector in 2006 between 
June and November has been searched for up-going neu- 
trino candidates. The search proceeds by a series of cut 
levels intended to remove down-going events as shown in 
Table HI Initially, hit cleaning is applied which removes 
all DOM hits which fall out of a 4 fxs time window, and 
all DOM hits without another DOM hit within a radius 
of 100 meters and within a time of 500 ns. After hit 
cleaning, we re-trigger, insisting that at least 8 DOM hits 
survive hit cleaning. Simple first-guess reconstruction al- 
gorithms running at the South Pole were used to filter 
out clearly down-going events. Events with fewer than 
11 DOMs hit were also filtered to meet bandwidth con- 
straints from the South Pole. The remaining events were 
transmitted to the data center in the Northern hemi- 
sphere via satellite and constitute the filter level of the 
analysis. At the data center, we reconstructed the direc- 
tion of events using a maximum-likelihood technique sim- 
ilar to the AMANDA muon reconstruction [!§]■ Events 
which were reconstructed as down-going were discarded. 
Despite the fact that remaining events appear up-going, 
the data is still dominated by misreconstructed down- 
going events. These down- going events are removed by 
additional quality cuts. Events which pass these quality 
cuts constitute the neutrino candidate dataset. 

Simulated events fall into the three categories shown in 
Table |TJ "Single shower" events arise from single cosmic- 
ray air showers in the atmosphere above IceCube and 
result in a single muon or bundle of collinear muons in 
IC-9. "Double shower" events come from two uncorre- 
cted air showers which happen to occur within the 5 /j,s 
event window. The CORSIKA[l3| air shower simulation 
program was used for the simulation of down-going sin- 
gle and double air-shower events. Finally, "atmospheric 
neutrino" events are muon neutrino events from pion and 
kaon decay in air showers in the northern hemisphere. 
The atmospheric neutrino model of [14j | and its exten- 
sion up to TeV energies fill ] as well as the cross-section 
parametrization of [161 ] were used to model the up-going 
muon rates due to atmospheric neutrinos. 

The events which are reconstructed as up-going are 
completely dominated by down-going muons from single 
and double-shower cosmic-ray events. Misreconstructed 
events are typically of low quality as measured by two pa- 
rameters, the number of direct hits Ndir, and the direct 
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FIG. 1: Data vs Cut Strength. Shown is the remaining num- 
ber of events as the cut strength S C ut (defined in the text) is 
varied. Curves are shown for the data and the total simulation 
prediction. Also shown is the prediction due to atmospheric 
neutrinos alone. The selection from the text corresponds to 
a cuts strength of S cu t = 10, and is denoted by an arrow. At 
this point, the data are dominated by atmospheric neutrinos. 



length Ldir- A direct hit is a photon arrival in a DOM 
which is detected between -15 and +75 ns of the time 
expected from the reconstructed muon with no scatter- 
ing. Ndir is the total number of direct hits in an event. 
The direct length Ldir represents the length of the recon- 
structed muon track along which direct hits are observed. 
An event with a large number of direct hits and a large di- 
rect length is a better quality event because the long lever 
arm of many unscattered photon arrivals increases con- 
fidence in the event reconstruction. The strength of the 
quality cuts can be represented by a dimensionless num- 
ber Scut which corresponds to cuts of Ndir > Scut and 
Ldir > 25 • S cu t meters. In addition to these quality cuts, 
we impose a cut requiring that events have no more than 
46 DOMs hit, which eliminates only about 1% of the final 
event sample. The purpose of this cut is to leave the high- 
multiplicity data blinded for an anticipated search for a 
high-energy diffuse neutrino flux. Figure [T] shows how 
many events remain as we turn the cut strength up and 
increase the signal-to-noise ratio. The accurate simula- 
tion of mis-reconstructed down-going events requires ex- 
cellent modeling of both the depth-dependent ice proper- 
ties and DOM sensitivity. In this initial study we observe 
a 60%-80% rate discrepancy for misreconstructed events 
up to a cut level of about S cu t = 8 or so. Nevertheless, 
over more than four orders of magnitude, the background 
simulation tracks the data, the number of wrongly recon- 
structed tracks is reduced, and for S cut > 10, the data 
behaves as expected for atmospheric neutrinos. From 
simulation, we expect neutrinos with energies between 
about 10 2 and 10 4 GeV, peaked at 1000 GeV, to survive 
the analysis cuts. 

In 137.4 days of livetime we expect 211 ± 76.1(syst.) ± 
14. 5(siai.) atmospheric neutrino events to survive at 
Scut = 10 and 234 events are measured. Above a zenith 
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TABLE I: Event Passing Rates (Hz). Shown are the event passing rates through different processing levels for the simulated 
event categories and for experimental data. The trigger level comprises the events triggering the detector after hit cleaning 
and re-triggering. The filter level comprises events which passed the online filtering conditions. Rates are also shown for events 
which reconstruct as up-going with and without the final quality cuts applied (see the text for cut definition). Note that 
the rates from air-shower events have been multiplied by 0.90 so that the simulation and data agree at trigger level. This is 
consistent with an approximately 20% uncertainty in the absolute cosmic-ray flux. For the final sample, statistical errors are 
given for the data and systematic errors are given for the atmospheric neutrino simulation. 

Criterion Satisfied Data Single Shower Double Shower Atmospheric Neutrinos 

Trigger Level 124.5 124.5 L5 6.6xl0" 4 

Filter Level 6.56 4.96 0.45 3.7xl0~ 4 

Up-going (S^t = 0) 0.80 0.49 0.21 3.3xl0" 4 

Up-going (S cut = 10) 1.97 ■ 10~ 5 ± 0.12 ■ 10~ 5 - - 1.77 ■ 10~ 5 ± 0.63 ■ 10" 5 

Up-going (Scut = 10 and 9 > 120) 1.19 • 10 -5 ± 0.10 ■ 10 -5 - - 1.42 • 10~ 5 ± 0.51 • 10~ 5 



of 120 degrees, where the background contamination is 
small, we measure 142 events with an expectation of 
169 ± 60.9(syst.) ± 13.0(stat) due to atmospheric neu- 
trinos. The principal systematic uncertainty in this at- 
mospheric neutrino expectation is due to the approxi- 
mately 30% theoretical uncertainty in the atmospheric 
flux normalization The other significant systematic 
error is due to uncertainties introduced by the modeling 
of light propagation and the detection efficiency of Ice- 
Cube DOMs. The uncertainty in the atmospheric neu- 
trino rate due to this modeling is estimated at 20% and is 
obtained in this initial study by examining changes in the 
neutrino passing rate when varying the cuts to account 
for the background simulation disagreement in Fig. [TJ 

Figure [2] shows the measured zenith distribution for 
the final event sample along with the atmospheric neu- 
trino prediction. The zenith angle distribution agrees 
well with atmospheric neutrino simulation for vertical 
events above about 120 degrees. The observed excess is 
believed to be residual contamination from down-going 
single and double cosmic-ray muons. This excess disap- 
pears if we tighten the cuts beyond S cu t = 10, suggesting 
that the recorded events at the horizon are of typically 
lower quality than expected from atmospheric neutrino 
simulation. This reinforces the belief we are seeing resid- 
ual background at the horizon. Above about S cu t = 12, 
with low statistics, the data at the horizon are consistent 
with a pure atmospheric neutrino signal. Figure [3] shows 
the azimuth distribution with the IC-9 geometry in the 
inset. The azimuth distribution has two strong peaks 
corresponding to the long horizontal axis of the IC-9 de- 
tector. The cut of 250 meters on event length constrains 
near-horizontal events that can be accepted along the 
short axis of IC-9 since the string spacing is 125 meters. 
We expect more uniform azimuthal acceptance in future 
seasons as the detector grows and becomes more sym- 
metric. 

We can characterize the response of the detector to 
neutrinos with an effective area A e g which is a function 
of neutrino energy E and neutrino zenith angle 9. The 
function A e g(E,0) is defined as the function which sat- 
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FIG. 2: Distribution of the reconstructed zenith angle 6 of the 
final event sample. A zenith of 90 degrees indicates a horizon- 
tal event, and a zenith of 180 degrees is a directly up-going 
event. The band shown for the atmospheric neutrino simu- 
lation includes the systematic errors from the text, and the 
error bars on the experimental data are statistical. Note that 
uncertainty due to the atmospheric neutrino flux is an uncer- 
tainty in normalization and is nearly independent of zenith 
angle. 
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R = dE dtt-$(E,e)-A e g(E,( 



(1) 



where <&(E,Q) is an arbitrary diffuse neutrino flux and 
R is the corresponding rate of events surviving analysis 
cuts. Figure 0] shows the effective area of IC-9 to neu- 
trinos with the event selection of S cu t = 10, both for 
neutrinos near the horizon and for nearly vertical neutri- 
nos. The effective area to neutrinos is much smaller than 
the geometrical area of the detector, due to the smallness 
of the neutrino cross-section. Above 10 5 GeV, the Earth 
starts to become opaque to neutrinos, and the highest 
energy up-going neutrinos can only be detected at the 
horizon. 
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FIG. 3: Distribution of the reconstructed azimuth angle <j> 
the final event sample. The band shown for the atmospheric 
neutrino simulation includes the systematic errors from the 
text, and the error bars on the experimental data are statis- 
tical. The inset shows the horizontal locations of the strings 
making up IC-9 relative to the center of the future array. 
Note that uncertainty due to the atmospheric neutrino flux is 
an uncertainty in normalization and is nearly independent of 
zenith angle. 



IV. CONCLUSIONS 

In 2006, IceCube was approximately 10% deployed and 
acquiring physics-quality data. Atmospheric neutrinos 
serve as an irreducible background to astrophysical neu- 
trino observation, as a guaranteed source of neutrinos for 
calibration and verification of the detector, and may be 
studied as a probe of hadronic interactions at energies 
inaccessible to terrestrial labs. In the first 137.4 days 
of livetime we have identified 234 neutrino candidates 
with the IC-9 detector. For events above 120 degrees, 
this neutrino sample is consistent with with expectations 
for a pure atmospheric neutrino sample. Selection of 
events was done within six months of the beginning of 
data acquisition, demonstrating the viability of the full 
data acquisition chain, from PMT waveform capture at 
the DOM with nanosecond timing, to event selection at 
the South Pole and transmission of that selected data via 
satellite to the North. During the 2006-2007 season, 13 
more strings were deployed, bringing the total number 
of strings for the Inlce detector to 22. The deployment 
of IceCube will continue during austral summers until 
2010-2011, while the integrated exposure of IceCube will 
reach a km 3 • year sometime in 2009. 




log 10 (E/GeV) 

FIG. 4: The effective area of IC-9 to a diffuse source of muon 
neutrinos. The event selection from the text was used, except 
that there was no multiplicity requirement imposed on the 
simulated events. The effective area to muon neutrinos and 
muon antineutrinos have been averaged to produce the figure. 
The effective area is shown as a function of neutrino energy, 
both for vertical and horizontal neutrino events. Horizontal 
events have a zenith of 90 degrees and vertical up-going events 
have a zenith of 180 degrees. 
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